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This  work  uses  a  method  for  the  stochastic  reconstruction  of  catalyst  layers  (CLs)  proposing  a  scaling 
method  to  determine  effective  transport  properties  in  proton  exchange  membrane  fuel  cell  (PEMFC). 
The  algorithm  that  generates  the  numerical  grid  makes  use  of  available  information  before  and  after 
manufacturing  the  CL.  The  structures  so  generated  are  characterized  statistically  by  two-point  correla¬ 
tion  functions  and  by  the  resultant  pore  size  distribution.  As  an  example  of  this  method,  the  continuity 
equation  for  charge  transport  is  solved  directly  on  the  three-dimensional  grid  of  finite  control  volumes 
(FCVs),  to  determine  effective  electrical  and  proton  conductivities  of  different  structures.  The  stochastic 
reconstruction  and  the  electrical  and  proton  conductivity  of  a  45  (Jim  side  size  cubic  sample  of  a  CL,  rep¬ 
resented  by  more  than  3.3  x  1012  FVCs  were  realized  in  a  much  shorter  time  compared  with  non-scaling 
methods. 

Variables  studied  in  an  example  of  CL  structure  were:  (i)  volume  fraction  of  dispersed  electrolyte,  (ii) 
total  CL  porosity  and  (iii)  pore  size  distribution.  Results  for  the  conduction  efficiency  for  this  example  are 
also  presented. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  the  near  future  the  world’s  energy  demand  will  have  to  be 
supplied  using  energy  sources  different  from  fossil  fuels.  Within 
such  scenario  hydrogen  fuel  could  be  playing  an  important  role  as 
an  energy  vector  [1].  On  the  other  hand,  due  to  its  high  efficiency 
and  the  potential  of  being  used  in  a  great  variety  of  applications, 
proton  exchange  membrane  fuel  cell  (PEMFC)  is  being  a  much  stud¬ 
ied  and  developed  technology  in  the  last  two  decades  [2,3]. 

Within  a  PEMFC,  a  catalyst  layer  (CL)  is  a  small  and  thin  region 
where  electrochemical  reactions  take  place  as  well  as  many  trans¬ 
port  processes,  which  to  a  large  extent  determine  its  general 
performance  [4,5].  The  solid  phase  of  this  component  is  formed  by 
an  electronic  conductor  which  also  serves  as  a  support  for  a  nanos- 
tructured  catalyst,  which  should  be  in  contact  with  a  dispersed  ionic 
conductor  (ionomer)  to  complete  the  reactions  [6].  This  CL  needs 
to  be  porous  to  allow  the  access  of  reactant  gases  entering  the  cell, 
and  to  let  water  formed  by  reactions  exit  the  cell  either  in  a  vapor 
or  in  a  liquid  phase  [7,8].  The  porous  medium  also  provides  a  large 
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active  area  available  for  larger  reaction  sites  and  the  resultant  larger 
current  [8]. 

In  a  PEMFC  the  catalyst  is  typically  platinum  (Pt),  the  support¬ 
ing  material  most  commonly  used  is  carbon  (C),  while  the  ionomer 
is  a  polymer  named  nation®.  These  materials  will  be  referred  to 
as  CL’s  primary  components.  In  a  PEMFC  CL  the  reactant  gases  are 
transported  from  the  border  between  CL  and  the  so-called  gas  dif¬ 
fusion  layer  (GDL)  up  to  the  reacting  catalytic  sites  through  pores 
formed  during  fabrication  between  the  primary  components.  While 
the  catalyst  has  the  role  of  promoting  the  proper  electrochemical 
reactions,  the  carbon  collects  and  conducts  the  produced  electrons 
and  the  ionomer  should  conduct  protons  generated  or  consumed 
by  the  proper  reaction. 

Experimental  observations  at  the  micrometer  scale  suggest  that 
the  CL  is  formed  by  agglomerates  in  a  porous  matrix  [6].  The  concept 
of  agglomerates  has  been  developed  some  years  ago  for  electro¬ 
chemical  systems  and  it  has  been  extensively  applied  to  PEMFC 
[9-12].  Nevertheless,  obtaining  information  at  the  nanoscale  is  very 
difficult  using  present  experimental  characterization  techniques, 
which  do  not  allow  for  the  capture  of  images  at  that  scale  without 
modifying  the  original  structure  [13,14].  This  limitation  has  lead  to 
the  formulation  of  different  theories.  The  most  popular  is  the  the¬ 
ory  that  considers  agglomerates  as  small  spheres  covered  with  a 
thin  layer  of  ionomer  and  an  inside  of  Pt/C  particles  [15-18]. 
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The  structure  and  composition  of  a  PEMFC  CL  define  this  com¬ 
ponent  as  a  random  heterogeneous  material.  This  heterogeneity 
arises  from  the  fact  that  it  is  constituted  by  different  phases,  being 
“phase”  an  identifiable  domain  with  its  own  particular  properties 
that  differentiate  it  from  the  rest  of  the  other  phases  in  the  CL  (i.e. 
voids,  other  solid  material,  gases  or  liquids). 

Proportionality  coefficients  for  mass,  energy  and  charge  trans¬ 
port  in  a  heterogeneous  material  are  significantly  affected  by  the 
properties  of  the  phases  from  which  this  is  composed,  by  the  vol¬ 
ume  fraction  composition,  and  by  the  structure  of  such  phases.  For 
this  reason,  an  effective  transport  coefficient  (ETC)  is  defined  for 
a  heterogeneous  material  as  a  proportionality  coefficient,  which 
characterizes  the  domain  of  the  material.  For  a  randomly  formed 
heterogeneous  material  with  M  phases,  a  general  ETC,  I<e  is: 

/Ce=/(/C1,/C2,...,/CM;01,02,...,0M;^)  (1) 

where  subscript  of  variables  indicates  a  respective  phase,  I<  is  the 
proportionality  constant  for  that  phase,  0  is  the  phase  volume  frac¬ 
tion  and  £2  is  a  microstructural  information  of  the  domain  [19]. 

In  literature  there  are  different  mathematical  relationships 
to  determine  ETCs  [20-22].  For  example,  Maxwell-Garnett  [21] 
developed  a  relationship  to  determine  optical  properties  for 
diluted  dispersions  of  uniform  spheres.  Later,  Bruggeman  extends 
Maxwell’s  model  to  systems  with  random  dispersions  of  spherical 
particles  with  an  extended  range  size  [22]. 

A  more  recent  and  powerful  technique  to  determine  effective 
properties  of  random  heterogeneous  materials,  is  the  so-called 
“stochastic  reconstruction”,  which  first  applications  is  attributed 
to  Quiblier  [23].  This  technique  is  based  on  the  computational  gen¬ 
eration  of  a  mesh  that  characterizes  the  real  microstructure  of 
the  heterogeneous  material,  mathematically  described  by  statis¬ 
tical  functions  referred  to  as  “correlation  functions”.  Torquato  [19] 
offers  a  methodology  to  characterize  microstructures,  as  well  as 
fundamental  theory  to  estimate  effective  properties.  More  recent 
works  offer  modifications  to  Torquato’s  methodology  to  optimize 
the  quality  in  the  stochastic  reconstruction  or  to  reduce  conver¬ 
gence  time  during  computing  [24-27]. 

The  right  value  of  ETC  is  indispensable  for  an  adequate  numer¬ 
ical  simulation  of  systems  containing  one  or  more  components 
with  characteristics  of  a  heterogeneous  media.  An  incorrect  ETC 
value  will  alter  the  results  and  consequently  it  will  provide  a  wrong 
interpretation  of  the  transport  phenomena. 

On  the  other  hand,  PEMFC  phenomena  simulations  have  been 
of  great  interest  in  recent  years  due  to  the  obvious  advantages  over 
experimental  techniques:  (1)  lower  research  costs,  (2)  it  is  a  pow¬ 
erful  tool  for  the  design  of  more  efficient  devices,  and  (3)  it  offers 
the  possibility  to  better  understand  phenomena  within  the  fuel  cell 
in  often  unreachable  places. 

When  the  simulation  objective  is  the  prediction  of  gas  flow  field 
performance,  it  is  common  to  assume  the  CL  as  a  one-phase  element 
[28-30].  Nevertheless,  to  detail  the  functioning  of  a  CL,  a  mathemat¬ 
ical  approach  has  been  used  based  on  the  assumption  that  a  CL  is 
a  heterogeneous  medium  formed  by  agglomerates  [11,31-33].  In 
this  case,  Bruggeman’s  approach  is  the  most  popular.  The  stochas¬ 
tic  reconstruction  has  been  successfully  used  in  the  simulation  of 
a  PEMFC  CL.  Wang  et  al.  [34]  proposed  a  direct  numerical  simula¬ 
tion  (DNS)  as  a  way  to  predict  the  CL’s  transport  phenomena  that 
take  place  at  the  micrometer  scale.  Unfortunately,  the  large  num¬ 
ber  of  control  volumes  required  to  detail  the  main  elements  that 
form  the  CL,  in  a  computing  domain  that  covers  the  whole  electrode, 
makes  DNS  an  expensive  computation-wise  method.  Recently,  Kim 
and  Pitsch  [17]  proposed  a  “sphere-based  simulated  annealing” 
method,  as  well  as  the  “lattice  Boltzmann”  method  to  reconstruct 
the  microstructure  of  the  CL  in  a  PEMFC  and  to  determine  the  effec¬ 
tive  diffusivity  of  the  structure,  respectively. 


In  this  work  a  theoretical  study  of  a  PEMFC  CL  by  the  stochas¬ 
tic  reconstruction  of  its  microstructure  and  a  scaling  technique  is 
presented.  The  scaling  technique  allows  a  significant  reduction  in 
computational  resources  needed  for  simulation.  Scales  are  defined 
based  on  different  observable  structures  at  different  order  of  mag¬ 
nitude  scales.  The  reconstruction  technique  employed  allows  the 
use  of  information  available  before  and  after  the  manufacture  of 
CLs.  The  so  reconstructed  structures  are  characterized  with  a  two- 
point  correlation  function  and  its  pore  size  distribution.  Results 
obtained  for  the  effective  conductivity  of  theoretical  samples  with 
different  nafion®  content  as  well  as  different  total  porosity  and  pore 
size  distribution  are  also  presented.  This  work  includes  the  follow¬ 
ing  steps:  (1)  microstructure  reconstruction  at  the  different  scale 
domains;  (2)  scaling  strategy;  (3)  statistical  characterization  of  the 
generated  structures;  and  (4)  analysis  of  ETCs  by  solving  the  charge 
transport  equation.  This  technique  can  also  be  applied  to  evaluate 
the  different  structural  theories. 


2.  Theory  and  calculation 

In  this  section  we  describe  the  microstructural  information  and 
the  perceptible  scales  of  the  CL.  The  mathematical  analysis  that 
conditions  the  grid  generation  (the  statistical  characterization  and 
the  ETC’s  determination)  of  a  subdivided  domain  by  n  samples  of 
inferior  length  scales  is  also  presented. 


2.1.  Structural  information  and  specific  scales  of  CL  in  PEMFC 


Although  the  nanometric  scale  information  of  CLs  is  limited, 
structural  and  compositional  information  does  exist  and  it  can  be 
obtained  before,  during  and  after  manufacturing  the  electrodes. 
Design  parameters  that  are  used  and  controlled  before  MEA’s  man¬ 
ufacturing  typically  include:  (1)  platinum  load  ( yPt )  in  mgPtcrn-2; 
(2)  ionomer  load  {fin)  which  is  a  relation  of  nafion®  weight  to  total 
electrode  weight  and  (3)  platinum  to  carbon  weight  ratio  ( 6Pt ), 
which  is  characteristic  of  the  catalytic  material  synthesis.  From 
such  information  and  from  other  physical  properties  of  its  compo¬ 
nents  (densities),  volume  fractions  of  phases  in  CL  can  be  calculated. 

On  the  other  hand,  the  solid  phase  volume  (Vs)  of  any  mixture 
is  equal  to  the  total  sum  of  individual  phase  mass  (m)  divided  by 
the  density  (p)  of  each  i  element  in  the  mixture,  Eq.  (2): 


(2) 


Using  design  parameters  ( yPt ,  0Pt )  in  Eq.  (2)  gives  the  volume 

occupied  by  each  primary  element  in  the  solid  phase  of  the  CL.  Eqs. 
(3),  (4)  and  (5)  can  be  used  to  calculate  platinum  (VPt),  carbon  (Vc) 
and  nafion  ( VN )  volume,  respectively: 


Vpt 


_ Vs _ 

[(Av(l  +  0pt)/0pt  Pn{]  -  Av)  +  (1/ Ppt)  +  (l/#pt  Pc)]ppt 


(3) 


f  Vpt  ppt 

V  Opt  pc 


) 


(4) 


f  Av(i  +  QptXVpt  ppt)\ 
\  -  MOpt  pn  ) 


(5) 


where  pPt,  pc  and  p n  are  platinum,  carbon  and  nafion  densities, 
respectively. 

On  the  other  hand,  Eq.  (6)  relates  the  electrode  total  porosity 
(&t)  with  the  solid  phase  volume  (Vs)  and  the  electrodes  total 
volume  (Vj) 


(6) 


Other  important  structural  information  that  can  be  obtained  before 
manufacturing  the  CL  is  the  Pt  and  the  C  particle  sizes.  Pt  is  typi- 


1250 


R.  Barbosa  et  al.  /  Journal  of  Power  Sources  196(2011)  1248-1257 


Fig.  1.  Domains  use  to  characterize  a  PEMFC  CL  at  different  scales:  (a)  “mesoporous”  where  the  blue  color  represents  the  pseudo-solid  phase  while  black  the  empty  phase; 
(b)  “agglomerate  and  micropores”,  the  blue  color  represents  the  agglomerates  while  black  the  empty  phase;  (c)  “inside  an  agglomerate”  here,  the  Pt  (~5  nm,  in  red),  carbon 
agglomerate  (~50nm,  in  blue)  and  ionomer  (in  light  blue)  can  be  distinguished.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  the  article.) 


cally  supported  on  C  during  catalyst  material  synthesis,  i.e.  before 
manufacture  of  CL,  which  also  conditions  the  final  electrode  struc¬ 
ture.  The  CL  manufacturing  method  determines  to  certain  degree 
its  final  structure,  but  such  effect  is  not  considered  in  this  work. 
After  manufacturing  the  CL  and  using  modern  techniques  such  as 
porosimetry  and  high  resolution  microscopy,  one  can  get  approx¬ 
imate  size  and  structure  of  agglomerates,  total  porosity  and  size 
pore  distribution.  In  this  work  we  assumed  that  physical  proper¬ 
ties  of  primary  components  are  not  modified  during  manufacture. 
This  assumption  is  valid  when  the  composite  material  is  fabricated 
using  physical  techniques  in  which  chemical  reactions  do  not  occur 
which  is  normally  the  case  in  state  of  the  art  MEA’s  manufacturing 
techniques.  Thin  layer  electrode  manufacturing  techniques  [35], 
are  among  the  most  popular  approaches  used.  Many  variants  to 
these  techniques  are  reported  in  literature  [36-40],  but  most  of 
them  involve  mixing  Pt/C  and  the  ionomer  in  a  liquid  solution  that 
is  physically  homogenized.  The  resultant  mixture  (catalytic  ink)  is 
then  deposited  on  a  substrate.  It  is  important  to  say  that,  regardless 
the  catalytic  ink  composition  and  type  of  substrate,  the  primary 
elements  are  effectively  distributed  randomly  (given  a  complete 
homogenization  stage). 

During  the  manufacturing  process  of  CLs  agglomerates  form 
composed  of  primary  components  plus  a  number  of  pores  of  a  vari¬ 
ety  of  sizes.  These  structural  features  of  the  CL  allow  us  to  define 
internal  substructures  at  different  scale  levels,  and  with  that  we 
can  establish  a  scaling  strategy  for  numerical  simulation. 

One  way  to  explain  these  apparent  structural  changes  at  differ¬ 
ent  scale  levels  is  by  picturing  yourself  at  infinite  getting  closer  and 
closer  into  the  CL  structure.  When  the  vision  scale  is  ~50,000nm, 
one  will  observe  a  homogeneous  dispersion  of  mesoporous  struc¬ 
ture;  for  a  scale  of  ~5000  nm  and  focusing  areas  where  mesoporous 
are  absent,  the  structure  is  defined  by  agglomerates  and  micro¬ 
pores;  finally  at  a  vision  scale  of  ~500nm  and  focusing  at  an 
agglomerate  the  observer  will  distinguish  a  structure  formed  by  a 
random  distribution  of  Pt/C  particles,  ionomer  and  probably  nano¬ 
metric  pores.  This  last  observation  actually  is  a  hypothesis  due  to 
the  fact  that  available  information  from  experimental  techniques 
is  limited  at  this  scale,  nevertheless  it  is  a  more  realistic  hypothesis 
compared  to  the  most  commonly  used  theory  where  agglomerates 
are  considered  spheres  with  Pt/C  particles  inside  and  covered  with 
a  thin  layer  of  ionomer  on  the  outside. 

In  Fig.  1  the  previously  described  scales  are  graphically  repre¬ 
sented  for  a  PEMFC  CL.  Fig.  1(a)  refers  to  the  so-called  “mesoporous” 


scale  with  a  ~50,000  nm  domain,  where  the  blue  color  represents 
the  pseudo-solid  phase  while  black  the  empty  phase;  Fig.  1(b)  rep¬ 
resents  the  “agglomerate  and  micropores”  scale  with  a  ~5000  nm 
domain,  here  the  blue  color  represents  the  agglomerates  while 
black  the  empty  phase;  finally  Fig.  1(c)  refers  to  the  “inside  an 
agglomerate”  scale  with  a  ~500  nm  domain  where  the  primary  ele¬ 
ments,  Pt  (~5  nm,  in  red),  carbon  agglomerate  (~50  nm,  in  blue)  and 
ionomer  (in  light  blue)  can  be  distinguished. 

In  this  work  the  sample  that  simulates  the  inside  of  agglomerate 
is  the  only  one  that  contains  the  primary  elements,  while  samples  of 
ascending  scale  contain  only  the  sample  of  the  immediately  smaller 
scale  and  its  corresponding  porosity.  We  define  the  relative  poros¬ 
ity  as  the  volume  fraction  occupied  by  the  empty  spaces  at  that 
scale,  ignoring  that  the  pseudo-solid  phase  (or  smaller  scale)  also 
contains  some  porosity.  The  volume  fractions  of  the  solid  phase  of 
each  sample  at  different  scales  can  be  related  mathematically  by 
Eq.  (7), 

n 

Sr  =  (7) 

1 

where  Sr  is  the  volume  fraction  of  the  solid  phase  in  the  entire 
CL,  Sj  is  the  individual  volume  fraction  of  the  solid  phase  in  each 
of  the  samples.  By  substitution  of  Eq.  (6)  in  Eq.  (7)  we  can  obtain 
a  mathematical  relation  that  calculates  the  “total”  porosity  in  the 
complete  CL  with  the  relative  porosities  of  each  of  the  samples  (Eq. 
(8)), 

n 

=  1  -  Pf1  -  (8) 

1 

where  represents  the  total  porosity  of  the  complete  CL  and 
the  relative  porosity  of  each  individual  samples. 

2.2.  Structural  reconstruction  algorithm 

Although  the  “annealing  method”  [19,24-27]  can  be  applied,  for 
simplicity,  our  reconstruction  method  considers  only  one  filter:  the 
one-point  correlation  function  (volume  fraction).  Fig.  2  presents  a 
general  flow  diagram  of  the  microstructural  reconstruction  algo¬ 
rithm.  When  the  structure  of  a  sample  is  composed  of  more  than 
two  phases,  the  subroutines  (dotted  line  square)  repeat  for  each 
one  of  them,  in  such  a  manner  that  the  structure  is  generated  in 
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Fig.  2.  General  flow  diagram  of  the  microstructural  reconstruction  algorithm.  The 
subroutines  in  the  dotted  line  square  are  repeated  for  each  reconstructed  phase. 


stages,  where  each  stage  corresponds  to  the  complete  generation 
of  a  phase. 

The  input  variables  in  the  microstructural  reconstruction  algo¬ 
rithm  are:  (i)  platinum  load  (ypt);  (ii)  nation®  load  {Pn)‘,  (iii)  weight 
percent  of  carbon-supported  platinum  ( 0Pt );  (iv)  total  porosity  and 
relative  porosity  of  n  -  1  samples;  (v)  domain  and  control  volume 
dimensions  at  each  scale;  and  (vi)  average  size  of  representative 
elements  at  each  scale. 

The  domain  is  fractioned  in  finite  control  volumes  (FCVs),  where 
each  volume  is  identified  by  an  index  number  randomly  distributed 
and  computer-generated  by  a  random  number  generator  [41  ]. 

The  reconstruction  method  of  each  phase  considers  three  forma¬ 
tion  stages:  (1)  the  center  of  the  representative  elements  at  each 
scale  (i.e.  primary  elements,  agglomerates,  mesoporous,  isles,  etc.) 
are  stochastically  distributed;  (2)  around  such  centers  and  in  one 
single  step,  a  previously  configured  specific  tridimensional  geom¬ 
etry  is  generated  (i.e.  amorphous  sphere  structures,  ellipses,  tubes, 
etc.);  (3)  in  a  random  manner  the  surroundings  are  filled  until  the 
required  volume  fraction  is  fulfilled. 

The  centers  of  such  structures  are  denominated  seeds  and  the 
number  of  seeds  distributed  in  the  computing  domain,  is  calculated 
using  Eq.  (9), 

NS  =  ^  (9) 

VLUi 

where  NS  in  the  integer  number  of  seed  control  volumes,  VCTl-  is  the 
number  of  total  control  volumes  necessary  to  fulfill  the  volume  frac¬ 
tion  of  the  representative  element  in  the  whole  computing  domain 
and  VCw  the  number  of  control  volumes  occupied  by  one  single 
representative  element.  The  index  function  (/s)  that  distributes  the 
seeds  has  the  following  general  form: 


f  1,  if  l<x<  NS 
1  0,  otherwise 


(10) 


The  structure  of  representative  elements  is  previously  defined  and 
designed  pixel  by  pixel  by  a  subroutine  according  to  the  known 
morphological  features  of  that  element,  to  then  be  systematically 


generated  around  all  seeds.  The  growth  of  the  surroundings  of  rep¬ 
resentative  elements  is  controlled  based  on  volume  fractions  as 
described  by  Eq.  (11): 


Wr-Va\ 

Vr 


(11) 


where  E  is  the  volume  error  that  exist  during  the  assignation,  VA  is 
the  assigned  volume  and  VR  the  reference  volume  for  the  specific 
phase.  VR  is  obtained  for  each  reconstructing  phase.  At  the  “inside 
an  agglomerate”  scale,  VR  is  obtained  by  Eqs.  (3),  (4)  and  (5)  for  plat¬ 
inum,  carbon  and  ionomer,  respectively.  At  the  “mesoporous”  and 
“agglomerate  and  micropores”  scales,  VR  is  obtained  by  Eq.  (6).  The 
range  of  available  adjacent  FCVs  (RID)  for  the  expansion  is  limited 
by  the  following  function, 

RID  =  VCT exp  (^2-)  (12) 

where  VCT  is  the  total  number  of  FCVs  in  the  sample  domain  and 
T  is  a  parameter  that  regulates  the  RID’s  value,  this  is  an  empiri¬ 
cal  relation  used  to  guarantee  that  VR  converges  to  VA.  The  index 
function  (Je)  for  this  assignment  is: 


Ie(NI)  = 


if  NJ  <  RID 
otherwise 


(13) 


where  NI  is  the  random  FCV’s  index  number.  After  the  expansion, 
the  structure  is  characterized.  The  two-point  correlation  function 
(S2)  of  a  homogeneous  medium  can  be  obtained  by  randomly  “toss¬ 
ing”  a  line  segment  of  r  length  with  a  specific  orientation  and 
counting  the  number  of  times  that  the  beginning  (x)  and  the  end 
(x  +  r)  of  the  line  fall  in  phase  j,  as  described  by  Eq.  (14): 

S2(x,  r)  =  (ij(x)Ij(x  +  r))  (14) 

where  the  angular  parenthesis  refers  to  the  statistical  average  when 
the  whole  domain  is  evaluated  and  Ij  is  equal  to  1  when  the  point 
belongs  to  phase;  as  described  by  the  following  identity  function: 


J  1,  if  x  is  in  the  phase  j 
1  0,  otherwise 


(15) 


Using  orthogonal  coordinates,  the  two-point  correlation  function 
employed  to  characterize  the  generated  structures  becomes: 


S2(r)  = 


1 

3N2 


N  N  N 

^S2il(r)  +  ^S2J(r)  +  J^S2<k(r) 

j,k=  1  k,i=  1  i,j=  1 


(16) 


where  N  is  the  computing  domain  cubic  length  and  S2iI(r)  is  the 
two-point  correlation  function  along  i  direction  (Eq.  (17)), 


N-r 

S2,i(r)  =  jv"T7X/('’J’  +  rj’ 

i=  1 


(17) 


Similarly  %  and  S2<u  can  be  defined. 

The  pore  size  is  determined  by  the  distribution  of  spheres  of 
different  radii  in  the  porous  structure,  where  the  pore  radius  is 
equal  to  radius  of  the  sphere  that  in  its  inside  is  formed  by  the 
“empty”  phase  [42].  To  avoid  that  large  spaces  are  fractioned  into 
smaller  spaces,  the  spheres’  radii  begin  to  be  modified  from  a  max¬ 
imum  limit  to  unity.  The  characterized  spaces  are  identified  by  the 
following  index  function, 


f  0,  if  B  c  X 

1  1 ,  otherwise 


(18) 


where  X  represents  the  porous  spaces  and  B  the  sphere.  The  maxi¬ 
mum  sphere’s  radius  is  determined  by  a  similar  function  evaluated 
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Table  1 

Dimensions  of  the  domain  and  control  volume  for  the  studied  scales. 


Sample 

Nx  =  Ny  =  Nz  (nm) 

dx-dy  =  dz  (nm) 

m3 

300 

3 

m2 

5000 

50 

Mi 

45,000 

300 

in  the  plane  i-j  along  k.  The  pore  size  distribution  in  each  sample 
at  different  scale,  is  defined  by  Eq.  (19), 


fra,Mn 


VPra,Mn 

VpMn 


(19) 


where  Fra?Mn  is  the  volume  fraction  that  characterizes  the  occu¬ 
pied  volume  by  pores  of  diameter  ra(VPrajMn )  over  the  total  volume 
occupied  by  the  empty  space  ( VPMn )  in  each  of  the  Mn  samples. 
For  scaling,  this  equation  can  be  extended  to  find  the  pore  size 
distribution  along  the  CL  (Eq.  (20)): 


The  effective  resistance  ( Reff )  of  a  heterogeneous  material,  com¬ 
posed  by  a  conductive  phase  and  one  or  more  insulating  phases, 
is  a  function  of  the  conductive  material’s  resistivity  (pm),  its  effec¬ 
tive  area  ( Ameff )  and  its  effective  length  ( Lmejf )  as  described  in  the 
following  equation, 

%  =  (23) 

Ameff 

From  a  different  angle,  the  Reff  can  be  calculated  as  a  function  of  the 
effective  resistivity  ( pmejf )  and  the  input  data:  area  of  the  sample 
(Am)  and  length  of  the  sample  (Lm),  as  described  by  Eq.  (24), 

%  =  /Vff ^  (24) 

Using  Eqs.  (23 )  and  (24)  we  can  obtain  pmeff  which  relates  the  effec¬ 
tive  resistivity  of  the  whole  sample  with  the  resistivity,  effective 
length  and  effective  area  of  conduction  of  the  conductive  phase, 
Eq.  (25), 


F7ya(Fra,Mn)  =  1  -  <*>,_,)  (20) 

where  FTra  is  the  volume  fraction  that  characterizes  the  occupied 
volume  by  pores  of  diameter  ra  over  the  total  occupied  volume  by 
the  empty  phase  in  the  whole  CL;  @T  is  the  CL’s  total  porosity  and 
the  sample’s  relative  porosity.  For  example,  when  n  =  1,  =  0. 

2.3.  Effective  ohmic  conductivity 

Ohm’s  law  relates  the  electric  current  (J)  of  a  conducting  mate¬ 
rial  directly  to  the  applied  voltage  or  potential  (AE),  as  described 
by  Eq.  (21), 

A  E 

J  =  kA ^  (21) 

where  k  is  the  proportionality  coefficient  so-called  conductivity 
which  is  a  property  of  the  material,  A  and  L  are  area  and  length  of 
charge  transport,  respectively.  By  applying  the  continuity  equation 
in  a  medium  discretized  by  FCV,  the  charge  conservation  transport 
equation  in  a  non-reactive  system  can  be  expressed  by  Eq.  (22): 

V  •  (/<mV0e)  =  0  (22) 

where  km  is  the  material’s  conductivity  and  0e  the  applied  poten¬ 
tial.  By  solving  the  charge  transport  continuity  equation  directly  in 
the  structural  mesh,  it  is  possible  to  obtain  the  average  of  all  local 
current  flows  (Ja),  generated  by  the  potential.  Taking  Eq.  (21)  and 
substituting  Ja,  we  can  obtain  the  effective  length  to  area  of  con¬ 
duction  ratio  (. LmefflAmeff=k((/)elJa ))•  As  it  will  be  shown  later,  this 
method,  will  simplify  the  structure  scaling. 


Table  2 

Porosity  of  simulated  structures. 


<M%) 

(%) 

(%) 

&m2  (%) 

0.50 

D1 

7.50 

45.95 

0.50 

D2 

10.00 

44.44 

0.50 

D3 

12.50 

42.86 

0.50 

D4 

15.00 

41.18 

0.40 

D1 

6.00 

36.17 

0.40 

D2 

8.00 

34.78 

0.40 

D3 

10.00 

33.33 

0.40 

D4 

12.00 

31.82 

0.30 

D1 

4.50 

26.70 

0.30 

D2 

6.00 

25.53 

0.30 

D3 

7.50 

24.32 

0.30 

D4 

9.00 

23.08 

0.20 

D1 

3.00 

17.53 

0.20 

D2 

4.00 

16.67 

0.20 

D3 

5.00 

15.79 

0.20 

D4 

6.00 

14.89 

Pmeff  — 


(25) 


As  the  effective  resistivity  characterizes  the  material  by  being  an 
intensive  property  and  being  able  to  be  extended  to  find  the  resis¬ 
tance  of  any  continuous  structure  formed  by  the  same  material,  this 
equation  can  be  generalized  to  find  the  resistivity  of  an  element 
formed  by  subdomains  of  smaller  scales  (Eq.  (26)): 


±±Ameffi±±Lmi 


i=l 


(26) 


where  pmeffT  is  the  effective  resistivity  of  the  global  domain  formed 
by  various  smaller  scale  subdomains  i;  Lmeff.  is  the  effective  length 
and  Ameff.  is  the  effective  area  of  the  phase  under  study  in  every 
subdomain.  Lmi  is  the  length  and  Am.  is  the  area  of  the  sample  of 
every  subdomain. 

To  normalize  and  generalize  results,  in  this  work  calculated 
resistivities  are  used  to  estimate  a  conduction  efficiency  (sk).  As 
conductivity  is  the  inverse  of  resistivity,  the  CL’s  effective  conduc¬ 
tivity  is  the  inverse  of  the  effective  resistivity  value  kmejfT  =  p“  . 
sk  is  calculated  by  comparing  the  effective  conductivity  with  the 
nominal  conductivity,  as  described  by  Eq.  (27).  By  substituting  Eq. 
(26)  in  Eq.  (27),  we  can  obtain  a  relation  that  provides  the  conduc¬ 
tion  efficiency  of  the  global  domain  ( ekj )  which  is  formed  by  several 
subdomains  of  smaller  scale  (Eq.  (28)), 


„  _  kmejf 
£k  ~  k 

Km 


(27) 


(28) 


For  example,  to  study  the  conduction  efficiency  with  three  scaling 
levels  the  equation  would  be: 

£j  —  f  Amefh Ameff 2Ameff3  \  /  ^m^m2^m3 
T  \  Lmeff i  Lmeff 2 Lmeff 3  )  \AmxAm2Am3 

Although  in  this  work  only  the  ohmic  conduction  efficiency  will  be 
determined,  other  transport  properties,  such  as  thermal  conduc¬ 
tivity  and  diffusion  coefficient,  just  to  mention  a  couple  of  them, 
could  be  determined  under  the  same  approach. 


3.  Simulation  conditions 


In  this  work  we  will  determine  the  electronic  and  ionic  conduc¬ 
tion  efficiency  of  a  CL  with  different  structural  features:  (1 )  ionomer 
electrode  load  (fiN)  in  the  range  of  20-80%wt;  (2)  CL  porosity  in 
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Fig.  3.  Two-point  correlation  function  (S2)  versus  a  non-dimensional  distance  (r/N) 
for  some  representative  structures  at  each  studied  scale,  (a)  shows  S2  for  the  carbon 
phase  at  the  M3  scale,  where  each  curve  corresponds  to  different  ionomer  loads 
(Av);  (b)  shows  S2  for  the  agglomerates  at  the  M2  scale;  and  (c)  shows  S2  for  the 
pseudo-solid  phase  at  the  M 1  scale.  In  (b)  and  (c)  each  curve  corresponds  to  different 
porosities  (<£,  ). 


Fig.  4.  Images  of  some  reconstructed  samples,  (a)  and  (b)  represent  the  M3  scale, 
with  fiN  equal  to  20  and  60%,  respectively,  the  blue  color  represents  the  carbon  phase 
and  dark  blue  the  ionomer  phase,  (c)  and  (d)  represent  the  M2  scale,  with  &m2  equal 
to  20  and  50%,  respectively,  the  blue  color  represents  the  agglomerates  and  dark 
blue  the  empty  phase,  (e)  and  (f)  show  the  Mi  scale  with  0Ml  equal  to  3  and  15%, 
respectively,  here  the  blue  color  represents  the  pseudo-solid  phase  and  dark  blue 
the  empty  phase.  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  the  article.) 


the  range  of  20-50%  and  (3)  pore  size  distribution  by  modifica¬ 
tion  of  the  relative  porosity.  All  generated  structures  have  the 
following  characteristics:  20%wt  Pt/C;  0.5 mg Pt cm-2;  average  Pt 
particle  diameter,  C  diameter  ( ERC ),  agglomerate  diameter  ( ERM ) 
and  mesopore  diameter  ( ERP )  of  3  nm,  50  nm,  500  nm  and  600  nm, 
respectively.  Density  of  materials  forming  the  solid  structure:  Pt, 
C  and  ionomer  are  21,450,  1800  and  2000  mg  cm-3,  respectively. 
Each  one  of  the  resulting  structures  is  generated  10  times  by  dif¬ 
ferent  random  series  and  then  their  values  are  averaged. 

Every  CL  is  then  studied  at  three  different  scales,  namely: 

(a)  Inside  an  agglomerate  (M3),  the  phase’s  reconstruction 
sequence  in  the  computing  domain  is:  (i)  carbon  agglomerate 
generation,  (ii)  platinum  on  carbon  placement-generation,  and 
(iii)  the  rest  is  ionomer. 

(b)  Agglomerates  and  micropores  (M2),  the  reconstruction 
sequence  is:  (i)  agglomerate  generation  and  (ii)  the  rest  is 
empty  room. 


(c)  Mesopores  ( M ! ),  reconstruction  sequence  is:  (i)  mesopore  gen¬ 
eration  and  (ii)  the  rest  is  a  pseudo-solid  structure. 

Table  1  briefs  the  dimension  of  the  computing  domain  (N*,  Ny 
and  Nz)  and  the  FCVs  (d*,  dy  and  dz )  for  each  of  the  scales  (M).  These 
values  are  based  on  average  diameters  of  representative  elements 
of  the  already  mentioned  scaling  method. 

In  this  work  the  pore  size  distribution  is  controlled  indirectly  by 
defining  the  relative  porosity  of  the  so-called  “mesoporous”  scale 
)  as  a  percentage  of  total  porosity  (Eq.  (30)), 


where  is  the  total  porosity  percentage  assigned  to  and 
is  the  total  electrode  porosity.  From  a  physical  point  of  view, 
00  is  directly  proportional  to  the  pore  size  distribution  because  it 
determines  the  volume  occupied  by  mesopores.  It  is  worth  men- 
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Fig.  5.  Pore  size  distribution  of  reconstructed  catalyst  layers,  (a),  (b),  (c)  and  (d)  correspond  to  20,  30,  40,  and  50%  of  total  porosity  (0t),  respectively.  The  arrows  show  how 
00  is  modified. 


tioning  that  the  relative  porosity  of  the  M2  scale  is  conditioned  by 
Eq.  (8).  Table  2  specifies  relative  porosities  used  in  this  work,  where 
Dl,  D2,  D3  and  D4  refer  to  equal  to  0.15,  0.20,  0.25  and  0.30, 
respectively. 

4.  Results  and  discussion 

The  present  methodology  proposes  the  study  of  the  following 
CLs:  4  porosities,  4  pore  size  distributions  and  13  different  ionomer 
loads  (from  20  to  80  in  steps  of  5%).  As  each  structure  will  be  gen¬ 
erated  by  ten  different  random  series  we  will  have  a  total  of  2080 
analyzed  structures.  Using  one  single  personal  computer  with  3  Gb 
of  RAM  memory  and  a  2.4  GHz  processor,  the  computing  time  is 
of  approximately  lOOh,  with  the  use  of  scaling,  not  including  the 
time  for  establishing  the  problem  and  for  the  analysis  of  results. 
The  scaling  method  not  only  reduces  computer  processing  time 
but  also  allows  a  detailed  study  of  the  CL  at  the  nanoscale  level. 
The  domain  lattice  of  the  cubic  simple  simulating  the  CL  )  is 
45,000  nm,  while  the  lattice  of  the  sample’s  control  volume  detail¬ 
ing  the  inside  of  an  agglomerate  (M3)  is  3  nm.  Without  the  scaling 
technique,  the  computing  domain  would  need  to  be  generated  by 
3.375  x  1012  control  volumes,  which  would  make  it  impossible  to 
do  with  the  previously  mentioned  computer. 

4.1.  Characterization  of  reconstructed  structures 

Fig.  3  shows  the  two-point  correlation  function  (S2 )  obtained  for 
some  representative  structures  at  each  studied  scale.  These  values 
are  plotted  versus  a  non-dimensional  distance  r/N,  where  r  is  the 
distance  between  two  points  and  N  the  domain  lattice.  In  Fig.  3(a) 
S2  is  shown  when  the  carbon  phase  is  analyzed  at  the  M3  scale.  Each 
curve  corresponds  to  a  different  ionomer  loading  (Pn).  In  Fig.  3(b) 
and  (c),  S2  is  calculated  for  the  “pseudo-solid”  phase  at  the  M2  and 
Mi  scales,  respectively.  In  these  last  two  figures,  each  curve  cor¬ 
responds  to  different  relative  porosities.  When  r  =  0,  S2  represents 
the  volume  fraction  of  the  characterized  phase.  All  correlation  func¬ 


tions  (S2 )  decay  exponentially  down  to  the  squared  volume  fraction, 
where  its  value  is  independent  of  r. 

Fig.  4  shows  representative  images  of  reconstructed  scales  cor¬ 
responding  to  S2  in  Fig.  3.  Fig.  4(a)  and  (b)  correspond  to  the  M3 
scale,  with  a  p^  of  20  and  60%wt,  respectively,  where  the  blue  color 
represents  the  carbon  phase  and  dark  blue  the  ionomer  phase. 
Fig.  4(c)  and  (d)  correspond  to  the  M2  scale,  with  @m2  of  20  and 
50%,  respectively,  where  the  blue  color  represents  the  agglomer¬ 
ates  and  dark  blue  the  empty  phase.  Fig.  4(e)  and  (f)  show  the  Mi 
scale  with  0Ml  3  and  15%,  respectively,  here  the  blue  color  rep¬ 
resents  the  pseudo-solid  phase  and  dark  blue  the  empty  phase.  In 
general  the  reconstructed  structures  look  alike  in  the  three  planes 
for  each  shown  condition.  Also  at  each  scale  it  is  clear  the  effect  of 
extreme  values  of:  ionomer  load  at  the  M3  scale  (Fig.  4(a)  and  (b)) 
and  relative  porosity  at  M2  (Fig.  4(c)  and  (d))  and  Mi  (Fig.  4(e)  and 
(f))  scales. 

The  pore  size  distribution  in  the  CL  is  obtained  using  Eq.  (20).  The 
resulting  expressions,  where  only  two  samples  contain  the  empty 
phase,  are: 

FTra(Fra,Ml)=%^(<Z>i)  (31) 

FTra(Fra,M2)=%^(02)(l-^l)  (32) 

Fig.  5  shows  the  pore  size  distribution  results  calculated.  Fig.  5(a), 
(b),  (c)  and  (d)  correspond  to  20,  30,  40,  and  50%  of  total  porosity, 
respectively,  for  a  hypothetical  CL.  In  these  plots  one  can  observe 
the  pore  size  variation  by  normalizing  From  that  figure  it  can 
be  seen  that  when  decreases,  the  frequency  of  the  ~50  nm  pore 
size  (left  peak)  increases  while  the  frequency  of  the  ~300  nm  pore 
size  (right  peak)  diminishes.  An  example  of  this  behavior  is  seen  in 
Fig.  5(a)  when  =  0.30  (D4),  FTr50  =  0.61  and  FTr300  =  0.28  and 
when  <£>0  =  0.15  (Dl),  F7Y50  =  0.71  and  FTr300  =  0.14.  The  arrows 
in  those  figures  point  it  out.  These  results  are  expected  as  0^  is 
directly  proportional  to  the  relative  porosity  of  sample  with  a  larger 
scale,  therefore  when  its  value  decreases  the  volume  of  those  pores 
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Fig.  6.  Schematics  of  the  general  boundary  conditions  for  the  charge  transport  sim¬ 
ulation  of  the  cases  studied  at  the  different  scales. 

with  larger  size  decrease.  These  results  confirm  that  the  pore  size 
distribution  of  the  CL  was  modified  by  the  control  of  Such 
distribution  shows  two  peaks  (at  0.05  p,m  and  0.3  pan),  which  are 
strongly  influenced  by  the  carbon  particle  size  and  the  mesopore 
size  both  introduced  as  input  data. 

4.2.  Electric  and  ionic  transport  coefficients 

Fig.  6  shows  the  general  boundary  conditions  (in  2D)  used  dur¬ 
ing  simulation  of  the  cases  studied  at  the  different  scales.  As  the 
structures  are  homogeneous,  the  side  of  the  cubic  domain  that  will 
be  assigned  to  the  “flow  current  input”  is  not  relevant,  but  clearly 
the  output  should  be  located  at  the  opposite  side.  At  the  input  side, 
a  1.1  V  potential  is  specified  exclusively  on  the  FCV  of  the  phase 
under  study,  while  in  the  faces  of  FCV  of  the  isolating  phase  the 
boundary  condition  is  a  zero  current  flow.  At  the  output  a  1.0  V 
potential  is  specified  only  for  the  phase  under  study,  therefore  hav¬ 
ing  a  0.1  V  gradient.  Other  sides  in  the  cubic  domain  are  specified 
with  a  zero  current  flow  condition.  It  must  be  mentioned  that  these 
limiting  values  have  the  only  purpose  of  determining  the  conduc¬ 
tion  efficiency  and  do  not  represent  actual  voltage  distribution  on 
the  different  scales  at  the  same  time. 

As  the  material’s  resistivity  does  not  affect  results  of  the  conduc¬ 
tion  efficiency  (Eq.  (28)),  the  phase  under  study  is  specified  with  an 
arbitrary  high  conductivity  (1000S  m-1 )  while  the  isolating  phase 
with  an  arbitrary  low  conductivity  (0.0001  S  m-1 ).  At  the  scale 
the  phase  under  study  is  the  pseudo-solid,  at  the  M2  scale  the  phase 
being  studied  is  the  agglomerates  and  at  both  Mi  and  M2  scales,  the 
insulating  phase  corresponds  always  to  the  empty  phase.  To  deter¬ 
mine  the  ionic  conduction  efficiency,  the  phase  under  study  at  the 
M3  scale  is  the  ionomer  and  the  insulating  phase  corresponds  to 
both  carbon  and  platinum.  To  determine  the  electronic  conduction 
efficiency,  the  phase  under  study  at  the  M3  scale  corresponds  to 
carbon  and  platinum,  while  the  insulating  phase  is  the  ionomer. 

Fig.  7  shows  the  potential  distribution  in  volts  within  the  elec¬ 
tronic  conduction  phase  for  reconstructed  CLs  shown  in  Fig.  4.  The 
potential  distribution  is  estimated  using  the  boundary  conditions 
previously  mentioned,  while  the  current  flows  are  used  exclusively 
to  calculate  the  conduction  efficiency.  The  domains  of  Fig.  7(a)-(f) 
are  limited  by  an  imaginary  line  and  the  potential  distribution  is 
plotted  exclusively  in  the  electronic  conduction  phase. 

Fig.  8  plots  proton  (a)  and  electronic  (b)  conduction  efficiency 
versus  different  ionomer  loadings  in  reconstructed  CLs  with  30% 
total  porosity,  where  each  curve  represents  a  different  pore  size 
distribution.  In  that  figure  it  can  be  observed  that  the  proton 


Fig.  7.  Potential  distribution  (V)  of  the  electronic  conduction  phase  of  some  recon¬ 
structed  samples,  (a)  and  (b)  show  potential  distribution  for  the  M3  scale,  with 
equal  to  20  and  60%,  respectively,  (c)  and  (d)  show  the  potential  distribution  for  the 
M2  scale,  with  &m2  equal  to  20  and  50%,  respectively,  (e)  and  (f)  show  the  potential 
distribution  for  the  Mi  scale  with  0Ml  equal  to  3  and  15%,  respectively. 


conduction  efficiency  increases,  while  the  electronic  conduction 
efficiency  diminishes,  both  in  an  exponential  fashion,  when 
increases.  This  is  an  expected  behavior  as  when  grows  what 
really  increases  is  the  volume  fraction  of  the  ionomer  phase  while 
the  electronic  conduction  phase  decreases  inversely.  On  the  other 
hand,  it  is  relevant  to  point  out  that  the  pore  size  distribution  does 
not  affect  the  effective  ionic  or  the  effective  electronic  conduction 
for  the  same  pN  value. 

In  Fig.  9  we  show  results  obtained  for  proton  (a)  and  elec¬ 
tronic  (b)  conduction  efficiency  versus  different  total  porosity  in 
reconstructed  CLs  now  with  35%  of  ionomer  loading.  Each  curve 
represents  a  different  pore  size  distribution  and  it  can  be  noted  that 
both  efficiencies  electronic  and  ionic  conduction  decrease  almost 
linearly  when  the  CL’s  porosity  increases.  This  trend  is  caused  by 
the  smaller  volume  proportion  of  solid  components  (electronic  and 
ionic)  when  the  porosity  increases,  which  clearly  affects  the  effec¬ 
tive  conductivity  of  the  heterogeneous  material.  Also  it  should  be 
noticed  that  the  slope  for  the  electronic  efficiency  is  different  from 
that  for  the  ionic  efficiency.  As  we  shall  see  later,  the  decay  rate 
depends  on  fiN.  It  can  be  noticed  from  these  figures  that  the  pore 
size  distribution  does  not  affect  the  effective  ionic  or  the  effective 
electronic  conduction  for  the  same  value. 

Finally,  in  Fig.  10  results  obtained  for  the  same  properties  as  a 
function  of  porosity  are  shown.  As  the  pore  size  distribution  does 
not  affect  these  results,  conduction  efficiency  dependent  on  the 
pore  size  distribution  were  averaged  to  be  included  in  this  analysis. 
Besides  the  before  mentioned  trends,  one  can  also  observe  that  the 
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Fig.  8.  Effect  of  ionomer  load  (/ % )  on  proton  conduction  (a)  and  electronic  con¬ 
duction  (b)  efficiency  of  the  reconstructed  CLs  with  30%  porosity  (@t)-  Each  curve 
represents  different  pore  size  distribution. 
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Fig.  10.  Effect  of  ionomer  load  (ft n)  on  proton  conduction  (a)  and  electronic  conduc¬ 
tion  (b)  efficiency  of  the  reconstructed  CLs.  Each  curve  represents  different  porosities 
(<£»r).  The  effect  of  the  pore  size  distribution  was  averaged. 


decay  rate  of  the  proton  conduction  efficiency  increases,  while  the 
electronic  conduction  efficiency  diminishes  when  increases. 

For  the  experimental  validation  of  the  present  method  the  fol¬ 
lowing  must  be  considered:  (i)  The  microstructure  reconstruction, 
at  each  scale,  should  be  made  by  a  procedure  that  statistically 


cpT/% 
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Fig.  9.  Effect  of  porosity  ( 0T )  on  proton  conduction  (a)  and  electronic  conduction 
(b)  efficiency  of  the  reconstructed  CLs  with  35%  ionomer  load  (ft N ).  Each  curve 
represents  different  pore  size  distribution. 


ensures  the  correct  representation  of  the  experimental  structure, 
(ii)  The  experimental  effective  transport  coefficients  should  be 
exclusive  to  the  CL.  This  information  is  not  currently  available  in  the 
literature.  Nevertheless  Boyer  et  al.  [43]  conducted  an  experimen¬ 
tal  study  where  the  proton  conductivity  of  a  pseudo-CL  is  similar 
to  the  results  presented  in  this  paper.  It  is  noteworthy  that  exper¬ 
imental  validation  is  a  current  work  in  our  laboratory  and  will  be 
published  in  the  near  future. 

5.  Conclusions 

A  technique  to  determine  the  effective  electronic  and  ionic 
conduction  of  a  catalyst  layer  in  a  PEM  fuel  cell  has  been  pro¬ 
posed.  The  scaling  approach  can  be  applied  to  determine  other 
effective  transport  parameters  of  heterogeneous  structures  in  CLs. 
The  mathematical  relationships  can  be  systematically  applied  to 
determine  effective  coefficients  of  other  processes  such  as  ther¬ 
mal  conduction  and  mass  diffusion.  The  proposed  mesh  generation 
algorithm  makes  use  of  available  information  before  and  after  the 
CL  fabrication  and  the  generated  structures  are  characterized  by 
two-point  correlation  functions  and  pore  size  distributions.  The 
scaling  approach  significantly  reduces  the  computing  hardware 
requirement  for  the  detailed  study  of  nanoscale  structures.  Finally, 
the  techniques  presented  in  this  work  can  be  applied  to  study 
already  existent  CL  microstructures  as  well  as  to  propose  new  and 
more  efficient  CLs. 
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